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Abstract 

The properties of the hollow beam produced by a conical lens are studied in detail. 
In particular, the impact of a rounded vertex is examined. It is shown that it could 
lead to drastic changes in the transverse distribution of the hollow beam, determined 
by the ratio between the transverse size of the incident beam and the size of the 
blunt area. An adequate choice for this ratio allows us to either minimize the losses 
or optimize the distribution symmetry. 
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1 Introduction 

Axicons are a family of cylindrical symmetrical optical systems that produce 
a line focus rather than a point focus from incident collimated beam[l]. There 
are several types of axicons, working either by reflection or by transmission, 
and being either converging or diverging, but the most common one is probably 
the conical lens. Such a lens is a cone made with a material of index n, and a 
basis perpendicular to its main z axis (Fig. 1). An incident collimated beam, 
propagating along the z axis, is deviated towards the main axis of the lens, 
with an angle P = {n — 1) a, where a is the base angle of the cone, a is usually 
a small angle, typically few degrees or less. Beyond the lens, two zones must be 
distinguished: in the first one, just after the lens, all deviated beams spatially 
coexist partially (in the hatched zone of Fig. 1), resulting in a diffraction- 
free beam. Beyond this zone, beams deviated along different directions do not 
overlap any more, and the intensity becomes distributed on a ring. Such hollow 
beams have recently regain interest, in the perspective of their use in optical 
trapping of atoms [2-4] and Bose- Einstein condensates [5]. 
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Fig. 1. Schematic behavior, in the geometrical approximation, of a colhmated beam 
through a conical lens C of base angle a. The propagation axis is the z axis, with 
its origin at the axicon vertex. After the conical lens, the beam is globally deviated 
towards the z axis with an angle /3. Point V is on the axicon vertex, while point X 
corresponds to the abscissa zq where a hole appears in the center of the beam in 
geometrical approximation. The hatched zone between V and X is the Bessel zone. 
For 2; > zq, the hollow beam has a radius tq and a thickness p. 

These applications require specific properties, as e.g. a zero intensity in the 
center of the beam, as predicted by the diffraction theory [6]. However, in most 
of the experimental realizations, the obtained pattern is more complicated: a 
usual observed defect is a drastic asymmetry of the distribution with a tail 
towards the center, possibly secondary rings and thus a non zero intensity 
at the center [2]. Practically, this leads to the need of a mask to clean the 
inner region of the hollow beam, and thus a loss of light power and additional 
diffraction effects, which are both major inconveniences in the experimental 
uses of such hollow beams. The source of these differences with respect to 
the theory is not clearly identified. As the conical lens vertex plays a central 
role in the transverse distribution of the hollow beam, it is natural to search 
for the origin of these defects in the properties of this point, in particular by 
taking into consideration its real shape. In the present paper, we study the 
modifications induced by the bluntness of the conical lens vertex on the hollow 
beam transverse shape. We show in particular that the ratio between the 
transverse size of the incident beam and the size of the blunt area determines 
the transverse distribution, and that an adequate choice for this ratio allows 
us to either minimize the losses or optimize the distribution symmetry. 



2 Geometrical propagation 



Fig. 1 describes schematically the geometrical beam propagation through a 
conical lens. The diffraction-free zone appears between the conical lens vertex 
V, and the point X where all deviated beams separate. The spatial distribu- 
tion of the light field in this zone is the superimposition of the interferences 
between deviated beams and diffraction on the conical lens vertex V. The 
transverse distribution of the resulting field follows a Bessel function [7]. Its 
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exact nature depends on the shape of the incident beam. For example, an 
appropriate incident Lagucrrc- Gaussian beam is able to generate a high-order 
Bcsscl beam of arbitrary order [8]. The main interest of Bessel beams is that 
they propagate without diffraction, i.e. their intensity transverse profile re- 
mains invariant upon propagation on macroscopic distances [9]. Bessel beams 
produced by axicons have been extensively studied, within the framework of 
numerous applications, such as second-harmonic generation [10], optical para- 
metric generator [11], and atom trapping [12]. 

Beyond point X, i.e. for z > zq, with zq — ri/P, where rj is the radius of the in- 
cident beam, all deviated beams separate (in the geometrical approximation), 
and a hole with radius Rq — p appears in the center of the beam. i?o = 1^^ is the 
distance between the z axis and the ray crossing the conical lens vertex, while 
p is the thickness of the ring (p — in the present case). It is also practical to 
define the radius tq of the ring, with tq = Rq—p/2 (Fig. 1). Such hollow beams 
are fundamentally different from hollow Lagucrrc- Gaussian modes. First, in 
the last ones, the circulation of the phase around the pattern center is equal to 
27r/, where / > is the angular index of the mode, while in the axicon hollow 
beam, it is equal to that of the incident beam, i.e. zero in the case of a plane 
wave or a TEMqo mode. This phase variation has no consequence on the inten- 
sity transverse distribution, but is critical in field dependent systems, e.g. if 
the hollow beam interferes with another beam. The second point concerns di- 
rectly the intensity distribution: the axicon hollow beam can be realized with 
arbitrary ratio between the radius and the thickness of the beam, as both 
quantities depend differently on z, while in Laguerre-Gaussian modes, this 
ratio is fixed by the Laguerre polynomial distribution and does not vary with 
z, except of course if it is transformed by an axicon[13]. Axicon hollow beams 
can be used in different applications, from laser machining[14] to optical traps 
for cold atoms[2-4] and Bose-Einstein condensates [5]. 



3 Gaussian beam propagation 

In practical realizations, incident beams are often produced by a laser, and 
thus are gaussian beams. In that case, the incident beam cannot be consid- 
ered as collimatcd, but is characterized by a waist wq where the beam spot 
size is minimum. To control the waist location and size, a convergent lens 
is usually inserted before the conical lens. In particular, it is easy with such 
a doublet to locate the waist after the conical lens. Fig. 2 shows the propa- 
gation of a coUimated incident beam through a convergent lens followed by 
a conical one, in the geometrical approximation. As for coUimated beams, 
an incident ray at the distance r from the z axis, is deviated with an angle 
(3, with /5 (r) = {n — 1) a + r/ f, where / is the focal length of the conver- 
gent lens. We can also define a distance Zq such that for 2; > Zq, a hole ap- 
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pears in the center of the beam. For thin lenses, a simple calculation leads 
to Zq = Ti {f — d) / (rj + (n — 1) af) where d is the distance between the two 
lenses. As in the single conical lens case, Rq and p, which now depend both on 
z, can be easily calculated. In particular Rq = {n — 1) az has the same value 
as in the single conical lens scheme. The point focus in F is transformed in a 
ring focus of radius Rq = {n — 1) a {f — d) , at location F' ~ F, where F is 
the focus point of the convergent lens (the interval between both points is in 
a^). If the gaussian properties of the incident beam are taken into consider- 
ation, all these quantities becomes indefinite, as the transverse expansion of 
a gaussian beam is theoretically infinite. However, for large incident beams 
(i.e. Wi » A, Wi being the incident beam size at 1/e^ and A the wavelength), 
the geometrical approach remains a good approximation. Therefore, the above 
expressions, where rj is replaced by Wi, can be used, keeping in mind that they 
concern 1/e^ intensity limits. 



Fig. 2. Schematic behavior, in the geometrical approximation, of a collimated beam 
through a converging lens L followed by a conical lens C. Point V is on the axicon 
vertex, while point X corresponds to the abscissa zq where a hole appears in the 
center of the beam. The hatched zone between V and X is the Bessel zone. 

This geometrical approach is unable to describe the thickness p of the ring. 
To evaluate it, it is necessary to take into account the diffraction of the beam 
during the propagation, and in particular on the conical lens vertex. Such a 
theory shows that when the convergent conical lens doublet is illuminated by a 
gaussian beam, the transverse distribution of the ring at the waist location, is 
also gaussian, symmetrical with respect to the r = Rq circle [6]. Its thickness 
at 1/e^ is pq = 1.65wo, comparable to the beam waist wq produced by the 
converging lens without axicon. The field amplitude of the wave after the lens 
is given by the usual Fresnel approximation to the Kirchhoff integral: 



where A is the wavelength and k = 2it/X is the wave vector. u{P) is the 
complex amplitude of the field in point P of the observation plane, Ui (M) the 
incident field in point M of the diffraction plane and r (M) the transmission 
function of point M. The integration is done on the whole diffracting surface 
5". ( is the optical distance between M and P, and thus e^'^^/C depicts the 
spherical wave expanding out from point M. The inclination factor I depicts 
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the angle of the incident and diffracted beams with the normal to the surface 
S. For small apertures and paraxial beams, it is approximated by / = — 2, so 
that Eq. 1 becomes: 



The transmission function depends of course on the considered optical ele- 
ments. It is tl (M) and tc (M) for respectively a lens and a conical lens: 



where r is the distance from point M to the optical axis, / is the focal length 
of the lens and n the optical index. Eqs 3 do not take into account a fixed 
phase shift induced by the thickness of the lenses. 

To compute the field amplitude after the conical lens and compare it with 
experimental measurements, we directly evaluate Eq. 2 through a standard 
integration software, without any additional approximation. Experimentally, 
we used a commercial conical lens from Bern optics, with a = 2° and a radius 
of 5 mm. It is in BK7 glass with an optical index n — 1.51. The incident 
gaussian beam, provided by a laser diode with A = 852 nm, has a waist 
Wi = 645 /im, so that it can be considered as collimatcd, and we have zq — 
Wi/j3 = 36 mm. The signal is detected through 2D and ID CCD cameras, and 
therefore we record the intensity of the pattern rather than its field amplitude. 
Thus the comparison between experimental and numerical results is performed 
through the intensity transverse distribution. 

Let us first ilhistratc the limits of the present model with a spectacular ex- 
ample: we consider the transverse distribution of the experimental beam de- 
scribed above, at a distance z — 115 mm = 3.2 zq after the conical lens. Fig. 3 
shows the experimental (dashed line) and theoretical (full line) intensity dis- 
tributions. While the ring diameters arc identical, it is clear that the present 
model is unable to reproduce the details of the distribution. In particular, the 
experimental profile shows a thin main ring, with inner contrasted secondary 
rings, looking like diffracting rings. On the contrary, the model predicts a 
unique thick ring, slightly modulated by diffraction. Thus it is clear that the 
diffracting elements are not treated correctly in the model. 






(3a) 



Tc (M) — exp {ik {n — l)r tan a) 



(3b) 
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Fig. 3. Experimental (dashed line) and theoretieal (full line) transverse intensity 
profile of the beam in z = 115 mm = 3. 2^0 after the axicon, obtained when the inci- 
dent beam, with a waist wi = 645 /itm, crosses a unique conical lens. The theoretical 
profile corresponds to the modelization of the axicon with a ponctual vertex. The 
vertical scale corresponds to the theoretical curve, when the total light power is 100 
mW. Units for the experimental curve are arbitrary. 

4 Modelling of a real cixicon 



To enhance the model, we take into consideration the bluntness of the conical 
lens vertex, by introducing an hyperbolic correction to the conical lens shape 
(fig. 4). We introduce the radius of curvature R of the hyperbole, so that the 
thickness e (r) of the lens is: 



e (r) — eo — Rtan^ {a) 



\^ ^ RHa.n' (a) 

where Cq is the thickness of the lens for R — 0, i.e. when the vertex is a 
point. To evaluate the radius of curvature -Rexp of our actual conical lens, we 
use a fit to adjust experimental and theoretical intensity distributions. The 
best fit, corresponding to Rexp = 3.5 mm, leads to a satisfying concordance 
between experiments and theory (Fig. 5). The remaining differences between 
both curves concern the amplitude of the secondary rings, but not their num- 
ber neither their location. We have also applied our model to non coUimated 
beams, i.e. when the conical lens is used with a convergent lens. As shown 
on Fig. 6, the transverse distribution of the beam after the doublet is well 
reproduced. Therefore the hyperbolic approximation of the vertex appears to 
be good enough to predict and optimize the use of an axicon to produce a 
hollow beam. 

Fundamental differences appear between the present configuration and the 
predictions of [6]. In particular, the beam size in the focal plane is almost 2.4 
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Fig. 4. Schematic representation of the axicon shape taking into account the flatness 
of its vertex. All parameters are those of Fig. 3, with i?exp = 3.5 mm. The origin of 
z is the conical lens vertex. 
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Fig. 5. Same as Fig. 3, except that the theoretical profile corresponds to the mod- 
elization of the axicon with the shape of Fig. 4. 

times larger than expected (Fig. 7a). If this broadening is effectively linked to 
the curvature of the conical lens vertex, we expect that it decreases when the 
incident beam waist increases. Fig. 7b shows that indeed, if Wi is increased 
significantly, the width of the ring becomes similar to that predicted in [6]. 

Fig. 8 shows how the ring thickness, measured at 1/e^ of its maximum inten- 
sity, evolves as a function of z. Typically, three zones can be defined, delimited 
by the values Z\ ~ 45 mm ~ 1.9 2:0 and Z2 — 75 mm ~ 3.1 2:0- Fot z < zi, the 
main peak of the distribution is in r = 0, as it is typical for the Bessel zone. 
For zi < z < Z2, the thickness is almost constant, with p = I.Gwq ~ po, where 
Wo is the minimum waist of the single convergent lens, in the focus plane, and 



7 



14-1 I 
12- 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
r(inrn) 



Fig. 6. Experimental (dashed line) and theoretical (full line) transverse intensity 
profile of the beam in z = 65 mm = 2.7zo after the axicon, obtained when the 
incident beam, with a waist Wi = 645 /xm, crosses a convergent conical lens doublet. 
The convergent lens has a focal length / = 100 mm, and is located at z = — 10 mm. 
The theoretical profile corresponds to the modelization of the axicon with a rounded 
vertex. Vertical scales as in Fig. 3. 

Pq is the value predicted in [6]. In the third zone, referred as expanding zone 
in the following, and corresponding to z > Z2, p grows almost linearly with z. 
This behavior is very different from that expected. In particular, the narrowest 
intensity distribution is not in the focus plane, but in the constant zone, just 
after the Bessel area. However, the smaller ring thickness corresponds to that 
predicted by [6] in the focus plane, although it occurs outside the focus plane. 

But the main novelties of our model concern the shape of the transverse distri- 
bution of the intensity. As shown in Fig. 9, very different shapes are obtained 
in the constant and expanding zone. In the constant zone, the external narrow 
ring is accompanied by several secondary rings located inside the main one, 
and separated by dark rings. In the immediate proximity of the Bessel zone, 
the minimum of the dark rings is almost zero (Fig. 9a), and it increases as the 
expanding zone is approach (Fig. 9b). At the edge of the constant and expand- 
ing zones (Fig. 9c), the contrast between secondary and dark rings vanish, so 
that in the expanding zone, the transverse distribution may be described by 
one single asymmetric ring with a tail towards its center. The resulting pat- 
tern in the expanding zone seems to be more clean, but it requires the use 
of a mask to suppress the intensity at the center of the ring, leading to new 
diffracting distortions of the pattern. On the contrary, in the constant zone, 
and in particular in the proximity of the Bessel zone (in 2^1), a mask adjusted 
to the first dark ring cleans the pattern to the single external ring, without ad- 
ditional distortion. The resulting pattern is a single narrow ring with two stiff 
sides. Thus the constant zone appears to be the most appropriate to produce 
clean hollow beams. 
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Fig. 7. The full lines show the theoretical transverse intensity profile of the beam 
in the focal plane of the convergent conical lens doublet for an incident waist 
Wi = 645 in(a) and Wi = 3 mm in (b). The dashed lines gives the gaussian 
distributions with a width po, predicted in [6]. Vertical scales as in Fig. 3. 

Another important aspect in the experimental realization of hollow beams is 
either the total light power inside the final beam, or the peak intensity of the 
ring. The first point depicts the losses introduced by the mask, while the second 
one is linked to the thickness of the ring. If the mask is placed at its optimal 
location zi, as described above, it introduces losses of 22%, corresponding to 
the part of power distributed in the secondary rings. This value have to be 
compared to that obtained when the mask is placed in the focus plane, as 
suggested in [6]. In this case, the losses depend naturally on the mask radius, 
and thus of the ring thickness. To obtain a thickness comparable to that in zi, 
we introduce 40% of losses. To reduce losses to 22%, it is necessary to increase 
the thickness by a factor 2.3: it is clear that the use of a mask in zi produces 
narrower rings with less losses. As a consequence, as it can be easily realized 
by comparing Fig. 7a and Fig. 9a, the peak intensity of the ring in almost 
four times larger in zi than in Zw So from the energetic point of view also, 
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Fig. 8. Intensity distribution thickness versus z. The width is measured at l/e^. The 
dashed hne is the value predicted in the focus plane in [6] . z^, is the abscissa of the 
focus plane. 

the technique described above leads to a substantial gain. 

In the whole preceding discussion, we have emphasized one given configura- 
tion, with a lens with focal length / = 100 mm, with a distance d = 10 mm 
between the lens and the axicon. The conclusions apply of course to other 
values of / and d. As an example, we illustrate in Fig. 10 the hollow beam 
obtained for a different set of parameters, using the above technique. The re- 
sulting pattern is a single ring with a stiff inner border, and losses introduced 
by the mask are measured to be 0.19±0.01, in good agreement with the theory. 



5 Conclusion 

We have shown that taking into account the curvature of the conical lens 
vertex leads to an accurate description of the beam transformation through a 
conical lens and a convergent conical lens doublet. Calculations are in excellent 
agreement with the experiments. We show that the use of a mask allows to 
obtain an excellent quality hollow beam, assuming that the mask is located 
adequately, just after the Bessel zone. The quahty of the hollow beam consists 
in a single regular ring, with an intensity going abruptly to zero on its inner 
side, and no residual intensity in the center of the ring. This enables the 
use of such beams in interferometric experiments, as e.g. in the realization 
of optical lattices[15]. In most of the applications, the present system will 
win in convenience if a second conical lens, located just after the Bessel zone 
of the first one, is added to the convergent conical lens doublet. Indeed, the 
second axicon will rectify the phase surfaces and fix the ring radius versus 
propagation. A third conical lens before the doublet may be added to control 
independently Rq and p[16]. 



10 




Fig. 9. Theoretical transverse intensity profile of the beam after the convergent 
conical lens doublet for the same parameters as in Fig. 6 for (a) z = 55 mm = 2.3zo, 
(b) z = 65mm = 2.7zo, (c) z = 75mm = S.Izq and (d) z = 95mm = 4.020. AH r 
axes have the same absolute scale, from r = mm to r = 2 mm. Vertical scales as 
in Fig. 3 



11 




Fig. 10. Experimental record of the intensity pattern obtained in zi for the converg- 
ing conical lens, with / = 500mm and d = 20mm. In (a), without the mask; in (b), 
with a 2 mm mask. 
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